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Abstract 

A recent chiral perturbation theory calculation of the in-medium quark condensate 
is extended to the isospin-asymmetric case of pure neutron matter. In contrast to 
the behavior in isospin-symmetric nuclear matter we find only small deviations from the 
linear density approximation. This feature originates primarily from the reduced weight 
factors (e.g. 1/6 for the dominant contributions) of the 27r-exchange mechanisms in pure 
| neutron matter. Our result suggests therefore that the tendencies for chiral symmetry 

restoration are actually favored in systems with large neutron excess (e.g. neutron stars). 
We also analyze the behavior of the density-dependent quark condensate (qq)(p n ) m the 
! chiral limit — > 0. 

PACS: 24.85. +p, 12.38. Bx, 12.39.Fe, 21.65.+f 

1 Introduction and framework 

The quark condensate |(0|gg|0}| is an order parameter of spontaneous chiral symmetry breaking 
in QCD. With increasing temperature the quark condensate decreases (or "melts"). For low 
temperatures T this effect can be systematically calculated in chiral perturbation theory. The 
\q ■ estimate T c ~ 190 MeV [1, 2] for the critical temperature, where chiral symmetry will eventually 
be restored, has been found. This extrapolated value is remarkably consistent with T c = 
(192 ± 8) MeV [3] obtained in recent lattice QCD simulations (modulo still persisting disputes 
1X2 '. between different lattice groups [4]). 

The chiral condensate | (qq) | drops also with increasing baryon density. Presently, it is 
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not feasible to study this phenomenon rigorously in lattice QCD due to the problems arising 
from the complex-valued fermion determinant at non-zero quark chemical potential. As an 
alternative, the density dependence of (qq)(p) can be extracted by exploiting the Feynman- 
Hellmann theorem applied to the chiral symmetry breaking quark mass term m q qq in the 
QCD Hamiltonian. The leading linear term in the nucleon density p is then readily derived 
by differentiating the energy density of a nucleonic Fermi gas, pM^ + 0(p 5 ^ 3 ), with respect 
to the light quark mass m q . This introduces the nucleon sigma-term = (N\m q qq\N) = 
m q dM^ /dm q = (45 ± 8) MeV [5] as the driving term for the density evolution of the chiral 
condensate. Following this leading linear density approximation one would estimate that chiral 
symmetry gets restored at (2.5 — 3)po, with p = 0.16 fm~ 3 the nuclear matter saturation 
density. 

Corrections beyond the linear density approximation arise from the nucleon-nucleon corre- 
lations which transform the nucleonic Fermi gas into a nuclear Fermi liquid. Because of the 
Goldstone boson nature of the pion, with its characteristic mass relation ~ m q , the pion- 
exchange dynamics in nuclear matter plays a particularly important role for the in-medium 
quark condensate. In a recent work [6] we have used in-medium chiral perturbation theory to 
calculate systematically the corrections to the linear density approximation. This calculation 
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has treated in detail the long- and medium-range correlations arising from l7r-exchange (with 
m^-dependent vertex corrections), iterated l7r-exchange, and irreducible 27r-exchange includ- 
ing also virtual A(1232)-isobar excitations, with Pauli-blo eking corrections up to three-loop 
order. It was furthermore necessary to estimate the quark mass dependence of a NN-contact 
term which encodes unresolved short-distance dynamics. Employing a recent computation of 
the NN-potential in lattice QCD [7] (at three different pion masses) we have found that the 
contact term has a negligible influence on the in-medium quark condensate (in agreement with 
ref.[8] which follows a somewhat different approach). As a result, we have obtained in ref. [6] a 
strong and non-linear dependence of the "dropping" condensate on the actual value of the pion 
mass m n . In the chiral limit = 0, chiral symmetry would seem to be restored already at 
about 1.5po- By contrast, for the physical pion mass = 135 MeV, the in-medium condensate 
stabilizes at about 60% of its vacuum value around that same density (see Fig. 2 in the present 
paper). 

Having found such pronounced deviations from the linear density approximation, with a 
hindered tendency towards chiral symmetry restoration in isospin-symmetric nuclear matter, it 
is a logical next step to explore the effects of the additional isospin degree of freedom present 
in nuclear many-body systems. Such a novel study is the purpose of the present short paper 
where we consider the extreme isospin-asymmetric case of pure neutron matter. Interestingly, 
we find that the same chiral pion-exchange dynamics which stabilizes the quark condensate in 
isospin-symmetric nuclear matter does not alter (in a significant way) its linear decrease in pure 
neutron matter (for densities p n < 0.35 fm 3 ). At the same time the behavior in the chiral limit 

— > changes prominently once a nonvanishing isospin-asymmetry is present in the nuclear 
medium. 

As in ref. [6] , our starting point is the Feynman-Hellmann theorem. It relates the in- medium 
quark condensate (qq)(p n ) 1° the quark mass derivative of the energy density of pure neutron 
matter. Using the Gell-Mann-Oakes-Renner relation m%f% = —m q {0\qq\0) one obtains for the 
ratio of the in-medium to the vacuum quark condensate: 



portional to = (N\m q qq\N) = m q dM^ /dm q comes from the non-interacting Fermi gas 
including the (relativistically improved) kinetic energy. We mention here that /„- denotes the 
pion decay constant in the chiral limit and m% stands for the leading linear term in the quark 
mass expansion of the squared pion mass. Interaction contributions beyond the linear density 
approximation are collected in the function: 



defined as the derivative of the interaction energy per particle E n {k n ) with respect to m\. Note 
that small (explicit) isospin breaking effects [9] (in interactions etc.) are not included in our 
calculation. The asymmetry in isospin is entirely given by the filled Fermi sea of neutrons (and 
the empty one for protons). 

2 Selected interaction contributions 

In this section we present analytical results for the contributions to the function D n {k n ) as given 
by various classes of In- and 27r-exchange diagrams, calculated up to three- loop order in the 
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energy density p n E n (k n ). At the level of individual diagrams pure neutron matter and isospin- 
symmetric nuclear matter differ only with respect to the (overall) isospin factor. In Table I 
several of these relative isospin factors are listed. Therein, we refer to the corresponding 
equation for D(kf) in ref. [6] and prescribe to substitute the Fermi momentum: kf — > k n . 
For completeness and clearness the remaining 27r-exchange contributions to D n (k n ) are better 
written out explicitely. We are following the enumeration scheme introduced in Sec. II of ref. [6] . 



Eq. in ref. [6] 
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isospin factor 


1/3 


1/3 


1/6 


-1/3 


1/6 


-1/3 


1/6 


1/6 


1/3 



Table I: Relative isospin factors for several In- and In -exchange contributions to D n {k n ). 

The irreducible 27r-exchange with only nucleons in the intermediate state leads to the fol- 
lowing contribution [10]: 
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with u = k n /m n and A the regularization scale. The three-body Fock term related to 2n- 
exchange with virtual A-excitation reads: 
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with the auxiliary functions Gs,t(x,u) written in Eqs.(7,8) of ref. [11]. The dominant two-body 
terms scaling reciprocally with the delta-nucleon mass splitting A = 293 MeV take the form: 
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For the remaining two-body terms with a more complicated A-dependence we employ the 
spectral-function representation [11] and differentiate directly the imaginary parts of the nNA- 
loop functions with respect to m 2 . This gives: 
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after subtracting a term linear in the neutron density p n = k^/3iT 2 . The associated subtraction 
constant includes also pion-loop contributions with a nonanalytical dependence on the quark 
mass m q . We reinstore these distinguished pieces by the term: 
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Finally, there is the 27r-exchange two-body term generated by the 7T7riViV-contact vertex propor- 
tional to the low-energy constant C\ = —a N /Am^ + 0{m. K ) (measuring explicit chiral symmetry 
breaking in the 7riV-interaction). Its contribution to the function D n {k n ) reads: 

~ „ ,tr- ?\ Zgl^mif, , 4N 27 + 49w 2 , . ?N 27 71m 93m 3 1 , . 

Dn(k n )^ = J^{(14« 2 + 3 M 4 )arctan M +^-^ln(l + M 2 )--- — -^}. (8) 

In comparison to Eq.(20) in ref.[6] only the numerical coefficient of the last -u 3 -term has changed. 
This comes from the different weighting of the Hartree and Fock contributions in neutron matter 
as compared to isospin-symmetric nuclear matter. 

The sum of all terms, Eqs. (3-8) together with those obtained via the relative isospin factors 
in Table I, comprise the m 2 -derivative of the complete set of in-medium Itt- and 27r-exchange 
processes up to three-loop order in the energy density (with inclusion of explicit A (1232) degrees 
of freedom). 



3 Results and discussion 

We are using consistently the same parameters in the chiral limit as in our previous work [6], 
namely: f w = 86.5 MeV, g A = 1.224, Cl = -0.93 GeV" 1 and M N = A = 882 MeV. For the 
quark mass dependence of the short- distance dynamics (not controlled by the underlying chiral 
effective field theory) we adopt the result derived in ref. [6] via the short range part (r < 0.6 fm) 
of the NN-potential from lattice QCD [7]. Its effect on the in- medium condensate (qq)(p n ) is 
then again negligibly small. 

Collecting all the pieces entering into Eq.(l), the condensate ratio (<?<?) (p n )/(0|gg|0) in pure 
neutron matter comes out as shown by the full line in Fig. 1. The dashed line therein corre- 
sponds to the linear density approximation using the empirical central value of the nucleon 
sigma-term, <Jn = 45 MeV [5]. In contrast to the behavior in isospin-symmetric nuclear matter 
(redisplayed in Fig. 2 for comparison) the chiral pion-exchange dynamics in pure neutron mat- 
ter generates only small deviations from the linear dropping of the in-medium condensate with 
density p n . This feature is straightforwardly explained by the reduced isospin weight factors of 
the 27r-exchange mechanisms (1/6 for the dominant contributions). Fig. 3 shows separately the 
effects of the five classes of interaction contributions for neutron densities p n < 0.35 fm~ 3 . They 
are consecutively added up in the sequence: linear density approximation — > l7r-exchange — > it- 
erated l7r-exchange — ► irreducible 27r-exchange — > 27r-exchange with virtual A(1232)-excitation 
— > chiral symmetry breaking Ci-term. One observes that the effects from 27r-exchange cancel 
here almost completely such that the total result (full line) lies close to the l7r-exchange ap- 
proximation (dashed-dotted line). This behavior is markedly different from the situation in 
isospin-symmetric nuclear matter (see Fig. 6 in ref. [6]). Note also that all the results discussed 
so far refer to the physical value of the pion mass, m w = 135 MeV. 

For the sake of completeness, we show in Fig. 4 the neutron matter equation of state resulting 
from our choice of parameters. Besides the In- and 27r-exchange contributions described in 
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Fig. 1: Ratio of the in-medium chiral condensate in pure neutron matter to its vacuum value 
at the physical pion mass, m n = 135 Me V. The dashed line corresponds to the linear density 
approximation using the empirical central value = 45 MeV [5]. 
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Fig. 2: Ratio of the in-medium chiral condensate in isospin-symmetric nuclear matter to its 
vacuum value. The dashed line shows the linear density approximation. 
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Fig. 3: Ratio between the in-medium chiral condensate in pure neutron matter and its vacuum 
value. The five classes of interaction contributions are consecutively added in the sequence: 
linear — > In — > iterated — > 2ir — > A — > c\. 



refs.[10, 11] and those proportional to Ci, it includes an adjusted short-distance term of the 
form: E n (k n ) ( - ad ^ = -1.04 GeV" 2 /c 3 - 18.4 GeV~ 4 k„. For not too high neutron densities, 
Pn < 0.2 fm -3 , our perturbative calculation reproduces fairly well the result of the Urbana 
group [12] based on a sophisticated many-body calculation. At higher neutron densities we get 
(as in ref. [11]) a stiffer neutron matter equation of state. This has to do with a repulsive k % n - 
term generated by the 27r-exchange three-neutron interaction. Note however that this /c^-term 
drops out when taking the derivative with respect to m 2 . 

Finally, we turn to the behavior of the density-dependent quark condensate (qq){p n ) in the 
chiral limit m n — > 0. As demonstrated in Sec. II F of ref. [6] the (singular) chiral logarithms 
ln(m n /X) from irreducible 27r-exchange and those from the m^-dependent vertex corrections to 
the l7r-exchange cancel exactly in the case of isospin-symmetric nuclear matter. This subtle 
balance does not work anymore for pure neutron matter. The following singular piece remains: 

D n (k n )\ m ^ = + 6^ - Hg 4 A ) In ^ . (9) 

In order to understand the physical origin behind this infrared singularity let us consider the 
simplified situation with the p-wave pion-nucleon coupling set to zero: qa = 0. In that case 
only the 27r-exchange interaction generated by the Weinberg- Tomozawa term (at second order) 
survives. It gives rise to an isovector (~ T\ • r 2 ) central NN-interaction with m^-derivative equal 
to: 

dW c (q) _ _L f 1 ^m^ y/4ml + q 2 q+ + q 2 

dm 2 ~ 647r 2 ^l2 n ^ q R 2^ 

The two-body Hartree term (left diagram in Fig. 5) in infinite neutron matter is proportional 
to its value at zero momentum transfer q = 0, thus exhibiting the chiral logarithm, —1/2 — 
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P n [fm 3 ] 

Fig. 4: Energy per particle E n (k n ) of pure neutron matter as a function of the neutron den- 
sity p n = k n /3n 2 . It includes an adjusted short- distance term of the form: E n {k n )^ ad ^ = 
— 1.04 GeV~ 2 k 3 — 18.4 GeV~^k\. The dashed-dotted line stems from the sophisticated many- 
body calculation of the Urbana group [12]. 






Fig. 5: Three-loop Hartree and Fock diagrams generated by the Weinberg-Tomozawa contact 
vertex at second order. Their combinatoric factor is 1/4. 



\n(m^/X). On the other hand the limit — > of Eq.(10) exists and it is proportional to 
1/2 — ln(g/A). Consequently, the two-body Fock term (right diagram in Fig. 5), involving the 
integral: dqq 2 (2k n - q) 2 (Ak n + q)[l/2 - ln(g/A)] = 4Jfe£[5 - 41n(2Jfe n /A)]/3, stays finite 
in the chiral limit. The coordinate-space "potential" associated with the logarithm ln(g/A) 
in momentum space actually behaves as r -3 (for r > 0). Its long-range character causes the 
interaction energy density for a homogeneous neutron-sphere of radius R to grow logarithmically 
with the system size: 
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In other words, the usual the thermodynamic limit N — > oo, R — > oo, with N/R 3 = constant, 
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does not exist for such a long-ranged interaction. In this context the quantity (t^) / f% plays the 
role of an isovector "polarizability" for the coupling and exchange of a pair of massless pions 
which causes the infrared singular behavior. Since the expectation value of r 3 does not vanish 
whenever the nuclear medium is asymmetric in isospin, the corresponding in-medium quark 
condensate will encounter an infrared singularity in the chiral limit m n — > 0. 

In passing we note that the result for the quark condensate in pure proton matter (ignoring 
the Coulomb interaction) is the same as that for pure neutron matter, shown in Figs. 1,3. 
Our results for the in-medium chiral condensate can also be interpreted in terms of a density- 
dependent effective nucleon sigma-term, crjv,eff(p) = CT jv,eff(p) + ( r 3) 2 CT yv,cfr(p)> with (r 3 ) = (Z — 
N)/(Z + N) the relative isospin-asymmetry. Its isoscalar piece <7^\ s (p) drops from the value 
(Tat ~ 45 MeV in vacuum to about 20 MeV at p = 0.35 fm~ 3 (see Fig. 7 in ref. [6]). The isovector 
piece cjveff(p)> which is maximally active in pure neutron or proton matter, brings it then 
back to almost the vacuum value in the density region p = 0.35 fm -3 considered here. As 
a consequence, the nuclear scalar mean-field associated in certain models [13, 14] with the 
in-medium condensate (qq)(p), would have an attractive isoscalar component as well as an 
attractive isovector component with approximately opposite density dependence. 

In summary, we have used in-medium chiral perturbation theory to calculate the density- 
dependent quark condensate (qq)(p n ) m pure neutron matter. We have found that the same 2n- 
exchange dynamics which stabilizes the chiral condensate in isospin-symmetric nuclear matter 
[6] does not alter (in a significant way) its linear decrease with the neutron density p n . This 
different behavior originates from the reduced weight factors of the 27r-exchange mechanisms in 
pure neutron matter. We can conclude that the tendencies for chiral symmetry restoration are 
actually favored in systems with large neutron excess (e.g. neutron stars). Possible observable 
consequences of this feature should be further investigated. 
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